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INTRODUCTION 

Seismic  refraction  surveys  are  a  valuable  tool  for  shallow  engineering 
site  characterization  purposes.  Standard  refraction  analysis  uses  first 
arrival  time  data  to  interpret  subsurface  layer  velocities  and  thicknesses. 

In  spite  of  their  usefulness,  refraction  surveys  suffer  from  several  well 
known  shortcomings  (for  example  blind  zones  and  velocity  inversions)  and  some 
not  so  well  known  difficulties,  such  as  thin  layer  effects.  In  this  study, 
synthetic  seismograms  were  computed  using  the  reflectivity  method  for  several 
thin  layer  models.  These  synthetic  record  sections  were  then  treated  as  field 
data  and  the  standard  interpretation  processes  were  performed.  The  inter¬ 
preted  subsurface  models  were  then  compared  with  the  true  geologic  model 
originally  input  to  the  reflectivity  code  to  determine  the  practical  problems 
associated  with  the  conduction  of  refraction  surveys  in  the  presence  of  thin 
layers. 

PREVIOUS  WORK 

Some  studies  of  the  thin  layer  problem  for  reflection  surveys  have  been 
done  in  recent  years  (Refs.  1,  2  and  3).  Studies  of  thin  layer  effects  have 
been  less  prevalent,  especially  for  shallow  refraction  surveys  used  in  engi¬ 
neering  site  characterization.  Several  papers  discuss  the  general  physical 
phenomenon  associated  with  thin  layer  propagation  as  important  to  seismic 
refraction  surveys  (Refs.  4-10).  References  4  and  5  present  the  results  of 
a  laboratory  seismic  model  study  of  refractions  from  layers  of  finite  thick¬ 
ness.  References  6  through  9  describe  the  results  of  calculational  studies 
of  thin  layer  headwave  propagation,  while  Reference  10  discusses  empirically 
derived  attenuation  coefficients  for  headwaves  propagating  in  thin  layers  as 
opposed  to  those  observed  for  propagation  in  thick  refractors  of  two  differ¬ 
ent  types  of  lithology. 

Before  further  discussion  of  thin  layer  refraction  effects,  a  defini¬ 
tion  of  what  is  meant  by  thin  is  required.  The  term  thin  is,  of  course, 
relative  to  the  wavelength  of  the  seismic  waves  propagating  within  a  layer. 
Reference  10  defines  thick  as  something  greater  than  a  few  wavelengths;  a 
thin  refractor  is  said  to  be  something  less  than  about  one-half  of  a 


wavelength.  As  Reference  10  points  out,  however,  these  magnitudes  are 
imprecise  and  will  also  depend  on  the  impedance  contrast  between  the  layer 
and  surrounding  material  and  the  observation  range. 

Refractions  may  become  unreliable  for  velocity  and  depth  determinations 
when  they  occur  with  wavelengths  which  are  large  compared  to  the  layer  thick¬ 
ness  (Refs.  4-9).  When  layers  become  thin  relative  to  wavelength,  an  inter¬ 
ference  condition  is  set  up  with  multiply-reflected  waves  within  the  layer 
interfering  with  the  headwave  propagating  within  the  layer.  This  results 
in  increased  attenuation  as  well  as  variations  in  phase  velocity  as  a  function 
of  wavelength  when  compared  to  headwaves  propagating  in  thick  layers.  This 
leads  to  various  phenomena  observed  in  refraction  recordings,  such  as  shing¬ 
ling  discussed  in  Reference  6. 


CALCULATIONS 

The  objective  of  this  study  was  not  so  much  a  reexamination  of  the 
physics  of  thin  layer  wave  propagation  but  a  determination  of  its  impact  upon 
the  results  of  shallow  refraction  surveys  performed  for  engineering  site 
characterization.  Reference  11  presents  a  field  example  of  a  shallow  refrac¬ 
tion  record  section  recorded  in  an  area  where  a  near-surface,  thin  high 
velocity  layer  overlies  a  fairly  thin  layer  of  lower  velocity  material,  which 
in  turn  overlies  a  thick  layer  with  approximately  the  same  velocity  as  the 
near-surface  layer.  This  subsurface  configuration  produced  a  refraction 
record  section  with  a  time  gap.  Consideration  of  geometric  ray  theory  alone 
would  predict  that  headwaves  from  the  lower  thick  high  velocity  layer  would 
never  be  observed  as  a  first  arrival.  In  reality,  however,  the  wave  propa¬ 
gating  through  the  near-surface  thin  layer  attenuated  rapidly  with  distance, 
particularly  for  the  lower  frequencies.  Eventually,  the  energy  traveling 
through  the  near-surface  thin  layer  dies  out  completely,  leaving  a  fairly 
wide  gap  in  time  before  headwaves  from  the  lower  thick  layer  arrive. 

In  the  alluvial  valleys  of  the  Western  United  States,  caliche  deposits 
often  form  near-surface  thin  layers  of  relatively  high  velocity.  Shallow 
refraction  surveys  are  often  performed  in  these  areas  to  obtain  material 
property  information  for  use  in  the  prediction  and  analysis  of  high-explosive 
test  data.  Often  record  sections  are  obtained  which  bear  a  strong  resemblance 
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to  those  presented  in  Reference  11.  A  typical  example  recorded  at  White 
Sands  Missile  Range,  New  Mexico,  is  presented  in  Figure  1.  If  thin  layer 
propagation  effects  are  responsible  for  the  time  gap,  as  shown  in  Figure  1, 
then  application  of  standard  interpretation  procedures  will  result  in  inac¬ 
curate  models. 

To  determine  the  influence  of  these  caliche  beds  upon  the  results  of 
shallow  refraction  surveys,  a  series  of  synthetic  refraction  record  sections 
were  computed  using  the  modified  reflectivity  method  (Ref.  12).  The  reflec¬ 
tivity  method  allows  the  computation  of  complete  seismograms,  including  sur¬ 
face  waves,  leaking  modes,  and  all  body  wave  phases.  The  method  has  been 
used  to  aid  in  the  interpretation  of  very  long  range  refraction  data  on  the 
scale  of  tens  to  hundreds  of  kilometers  (Refs.  13  and  14).  To  the  best  of 
our  knowledge,  the  method  has  not  been  used  before  to  aid  in  the  interpre¬ 
tation  of  engineering  scale  refraction  surveys. 

The  reflectivity  technique  was  first  developed  by  Fuchs  as  discussed  in 
Reference  15.  The  technique  uses  the  Thomson-Haskell  matrix  formalism  and 
numerical  integration  in  the  frequency  wavenumber  domain  to  solve  the  equa¬ 
tions  of  motion  for  a  layered  medium.  Kind,  as  described  in  Reference  12, 
later  extended  the  technique  to  compute  synthetic  seismograms  for  the  case 
of  a  point  source  buried  in  a  layered  medium.  The  introduction  of  attenua¬ 
tion  to  the  technique  shifts  the  poles  of  the  Rayleigh  waves  away  from  the 
real  axis  of  the  wavenumber  plane  allowing  numerical  integration  along  the 
real  axis,  thus,  enabling  the  computation  of  complete  synthetic  seismograms 
rather  than  body  waves  alone. 

The  AFWL  version  of  the  reflectivity  code  first  calculates  the  Green's 
function  or  impulse  response  for  a  layered  medium.  The  Green's  functions 
are  then  convolved  with  a  source  time  function  and  instrument  response  func¬ 
tion  to  produce  a  final  synthetic  seismogram. 

Synthetic  seismograms  were  generated  for  the  two  geologic  models  shown 
in  Figure  2.  These  models  are  idealized  approximations  to  the  geology  pres¬ 
ent  at  an  alluvial  high-explosive  test  site  in  New  Mexico.  The  thin  layer 
was  intended  to  simulate  a  caliche  bed  possessing  a  moderate  degree  of 
cementation.  The  velocities  shown  in  Figure  2  were  derived  from  shallow 
refraction,  uphole,  downhole,  and  crosshole  surveys.  To  determine  the 
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Figure  1.  Example  of  field  record  section  containing  time  gap  (recorded 
in  the  northern  Tularosa  Basin  on  White  Sands  Missile  Range, 
Mew  Mexico) . 
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effect  of  variations  in  thickness  of  the  caliche  layer  upon  the  synthetic 
record  sections,  the  thickness  of  the  caliche  layer  was  varied.  The  first 
set  of  calculations  was  performed  for  Model  A,  shown  in  Figure  2;  the  second 
set  for  Model  B.  Synthetic  record  sections  were  computed  for  Model  B  with 
thin  layer  thicknesses  of  0.0,  0.06,  0.12,  0.25,  0.50,  and  1  m.  For  the 
case  where  the  subsurface  half-space  was  included,  Model  A,  synthetic  record 
sections  were  computed  for  thin  layer  thicknesses  of  0.0,  0.06,  0.12,  0.25, 
0.50,  1,  and  2.35  m.  Synthetic  seismograms  were  computed  at  ranges  of  1.5, 
3,  4.5,  6,  7.5,  and  9  m. 


A  pure  explosive  source  was  arbitrarily  placed  at  a  depth  of  1  m  since 
the  code  does  not  allow  the  use  of  a  surface  source.  After  the  bulk  of  the 
computations  was  completed,  a  few  runs  were  repeated  wit1''  a  source  depth  of 
0.5  m.  No  qualitative  changes  in  the  observed  thin  layer  effects,  due  to 
the  change  in  source  depth,  were  obvious.  The  time  function  for  the  source 
consisted  of  a  step  function  with  a  second  order  polynomial  rise  time  of 
5  ms.  The  instrument  response  was  constructed  to  closely  approximate  the 
response  of  the  typical  shallow  refraction  seismic  recording  system.  The 
geophone  response  simulated  4.5  Hz  natural  frequency  units  damped  at  0.707 
critical.  The  recording  system  response  rolled  off  at  400  Hz.  A  total  of 
256  samples  was  computed  for  each  synthetic  seismogram.  The  sample  interval 
was  0.3  ms.  Some  runs  were  made  with  a  source  rise  time  of  2  ms.  The  reflec¬ 
tivity  code  includes  attenuation  effects.  For  these  runs,  all  P  wave  Q 
factors  were  set  equal  to  100,  S  wave  Q  factors  were  set  equal  to  44. 

An  estimation  of  the  wavelength  for  the  energy  propagating  in  the  thin 
bed  may  be  made  by  assuming  the  period  to  be  equal  to  the  rise  time  (5  ms) 
of  the  source.  This  yields  a  P  wavelength  of  3.5  m  for  energy  in  the  caliche 
layer  and  1.8  m  in  the  surrounding  alluvium.  For  the  high  frequency  source 
(2  ms),  we  have  a  caliche  wavelength  of  1.4  m  and  0.73  m  in  the  surrounding 
material . 


RESULTS  OF  THE  CALCULATIONS 

Vertical  and  radial  synthetic  record  sections  are  shown  in  Figures  3 
through  14  for  thin  layer  thicknesses  of  0.0,  0.06,  0.12,  0.25,  0.50,  and 
1  m  for  the  subsurface  model  (Model  B)  without  the  lower  half-space.  The 
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Figure  7 


Vertical  synthetic  record  section  for  model  B 
with  0.50  m  thin  layer. 
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arrows  in  these  figures  indicate  picks  of  the  first  arrival  time.  The 
source  rise  time  used  for  the  calculations  shown  in  these  record  sections 
was  5  ms. 

Examination  of  the  record  sections  in  Figures  3  through  14  reveals  the 
effects  of  changing  thin  layer  thicknesses  upon  the  appearance  of  the  syn¬ 
thetic  seismic  sections.  The  most  striking  effect  is  the  weakening  of  the 
first  arrival  in  the  waveforms,  particularly  on  the  vertical  component. 

After  the  thin  layer  reaches  0.12  m  thickness,  the  first  break  becomes  emer¬ 
gent  and  difficult  to  pick. 

Figures  15  and  16  illustrate  the  changes  in  the  waveforms  which  occur 
at  6  m  range  with  changes  in  the  thickness  of  the  caliche  layer.  A  well- 
defined  clear  first  arrival  is  apparent  on  the  0.0  m  thickness  thin  layer 
waveforms.  This  arrival  becomes  less  defined  as  the  thin  layer  thickness 
increases.  The  first  breaks  appear  to  be  stronger  on  the  radial  component. 

Some  synthetic  record  sections  have  a  qualitative  similarity  to  the 
field  record  section  in  Figure  1,  although  the  purpose  of  the  calculational 
effort  was  not  to  directly  model  the  Figure  1  records.  As  in  the  synthetic 
record  sections,  a  region  appears  in  the  field  section  (just  past  halfway 
through  the  spread)  where  the  first  arrival  becomes  very  emergent  and  then 
dissipates  completely.  In  the  field  section,  a  strong  arrival  resumes  at  a 
later  time  creating  the  characteristic  time  gap  in  the  travel  time  curves. 
The  synthetic  sections  do  not  extend  far  enough  in  range  to  show  clearly 
the  resumption  of  strong  first  arrivals  after  the  region  where  the  initial 
first  arrival  decays. 

All  subsurface  models  for  which  the  calculations  were  performed  are, 
of  course,  examples  of  velocity  inversions.  The  thin  layer  is,  in  each  case 
higher  in  velocity  than  the  underlying  material.  Purely  geometrical  ray 
theory  would  predict  then  that  the  material  beneath  the  thin  layer  would 
never  be  detected  by  stan  ’rd  refraction  interpretation  techniques.  The 
purpose  of  this  study  was  not  an  attempt  to  develop  a  technique  for  solution 
of  the  velocity  inversion  problem,  but  to  examine  the  departure  from  geo¬ 
metric  ray  theory  under  conditions  when  the  thickness  of  a  geologic  bed  is 
small  when  compared  to  the  wavelength  of  the  seismic  energy  used  in  a 
refraction  survey. 
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Figure  16.  Comparison  of  radial  synthetic  seismograms 
for  the  thin  layer  thickness  indicated, 
model  B,  6m  range. 
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To  gain  some  insight  into  thin  layer  effects  upon  the  results  of  shallow 
refraction  surveys,  the  synthetic  record  sections  were  treated  as  real  field 
data  and  first  break  times  picked  as  indicated  by  the  arrows  in  Figures  3 
through  16.  Travel  times  from  all  the  record  sections  are  shown  in  Figures 
17  through  23.  The  solid  lines  indicate  true  travel  times  computed  for  the 
models  for  waves  propagating  in  the  near  surface  material  and  the  thin  fast 
layer  assuming  geometric  ray  theory.  The  near  source  curvature  of  the  366  m/s 
curve,  as  well  as  the  non-zero  intercept  of  both  curves,  results  from  the 
1  m  depth  of  the  source.  The  data  points  are  the  first  break  times  deter¬ 
mined  from  all  the  synthetic  record  sections. 

Comparison  of  the  theoretically  predicted  travel  times  in  Figures  17 
through  23,  with  the  first  break  times  picked  from  the  synthetic  record  sec¬ 
tions,  reveals  large  discrepancies  between  the  observed  and  predicted  values 
for  the  thinner  layers.  There  is  a  good  correlation  between  the  observed  and 
geometrically  predicted  arrival  times  for  the  2.35  m  thickness  case.  The 
2.35  m  case  can,  however,  hardly  be  considered  a  thin  layer  case  since  the 
caliche  layer  lies  directly  on  top  of  the  lower  half-space  which  has  veloci¬ 
ties  near  those  of  the  caliche  layer. 

How  do  the  discrepancies  between  the  first  break  picks  and  the  true 
(geometrically  predicted)  travel  times  affect  what  the  interpreted  subsur¬ 
face  structure  will  be  compared  to  the  true  (that  which  was  input  to  the 
reflectivity  code)  structure?  To  determine  the  errors  which  would  be  present 
in  the  interpreted  subsurface  model  had  the  synthetic  record  stations  been 
real  field  data,  some  of  the  picked  first  break  times  were  input  to  an  auto¬ 
mated  seismic  refraction  interpretation  code.  This  code,  SIPT,  is  a  24- 
geophone  version  of  the  code  described  in  References  16  and  17.  The  code 
accepts  time  distance  data  from  refraction  surveys  and  uses  a  delay  time 
approximation  to  produce  a  first  cut  subsurface  model.  This  model  is  then 
adjusted  by  up  to  three  passes  of  ray  tracing  to  reduce  the  discrepancy 
between  observed  and  computed  arrival  times  at  each  geophone.  Output  of 
the  code  is  a  velocity-depth  structure  containing  up  to  five  layers.  The 
SIPT  code  is  used  as  part  of  the  standard  refraction  interpretation  process 
at  AFWl. 


22 


Comparison  of  true  versus  picked  arrival  Limes  for  no  thin  layer  cases 
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figure  19.  npp.iri:  i  i  true  versus  picked  arrival  times  for  0.12  m  thin  layer  cases 
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re  22.  Comparison  of  true  versus  picked  arrival  times  for  1  m  thin  layer  cases. 


23.  Comparison  of  true  versus  picked  arrival  times  for  2.35 


To  interpret  the  synthetic  record  section  arrival  times,  they  were 
treated  as  fully  reversed  times  (times  for  the  reverse  were  made  identical 
to  the  forward  times)  from  a  6-geophone  refraction  spread  with  a  geophone 
spacing  of  1.5  m.  The  code  requires  that  a  layer  number  be  assigned  to  each 
arrival  time.  For  each  case,  the  arrival  at  the  first  geophone  was  said  to 
be  associated  with  Layer  1;  all  other  arrivals  were  assigned  to  Layer  2. 

Results  from  the  interpretation  exercise  are  presented  in  Figures  24 
through  29,  starting  with  the  results  for  the  0.06  m  layer  in  Figure  24  and 
ending  with  the  2.35  m  layer  case  in  Figure  29.  Intuition  would  suggest 
that  a  layer  too  thin  to  be  detected  as  a  discrete  unit  would  simply  be 
averaged  in  some  thickness  weighted  fashion  with  the  velocities  of  the  sur¬ 
rounding  material.  In  some  of  the  interpretations ,  this  appears  indeed  to 
have  happened  (for  example,  the  0.06  m  layer  case  in  Figure  24,  perhaps 
Figure  25  which  is  the  0.12  m  case,  and  the  0.25  m  layer  case  in  Figure  26). 
For  the  thicker  thin  layer  cases,  0.50  m  in  Figure  27  and  1  m  in  Figure  28, 
the  interpreted  subsurface  models  are  clearly  not  anywhere  near  being  a 
thickness  weighted  average  of  the  thin  layer  and  surrounding  material.  For 
the  0.50  m  layer  model,  we  obtain  an  interpreted  section  with  a  velocity 
interface  at  approximately  2.3  m  depth,  a  significant  departure  from  the 
true  depth  of  1  m.  In  the  interpretation  of  the  1  m  thick  layer  shown  in 
Figure  28,  the  interpreted  velocities  are  near  the  true  values;  however,  the 
depth  of  the  velocity  interface  is  significantly  greater  than  the  true  value, 
and  the  code  has  placed  topography  on  the  interface  in  an  effort  to  fit  the 
first  break  times.  The  2.35  m  interpretation  shown  in  Figure  29  is  quite 
close  to  the  true  subsurface  model.  The  velocity  interface  is  at  the  right 
depth,  1  m,  while  the  velocities  are  near  the  true  values.  The  2.35  m  layer 
case  is  probably  not  properly  considered  as  a  thin  layer  case  since  the  2.35  m 
caliche  layer  is  in  direct  contact  with  the  underlying  half-space  with  almost 
the  same  velocities.  The  good  correlation  between  the  interpreted  and  true 
structures  for  the  2.35  m  case  verifies  the  interpretational  capabilities 
for  the  SIPT  computer  code. 
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of  the  synthetic  record  sections.)  The  travel  times  used  were  picked 
the  model  B  vertical  record  section. 
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DISCUSSION 


Interpretations  of  the  synthetic  record  sections  have  shown  that  the 
thin  layer  effects  can,  in  certain  case,  have  fairly  deleterious  effects  upon 
the  results  of  shallow  refraction  surveys.  Intuition  would  suggest  that  in 
cases  where  the  thin  layer  is  not  directly  detectable,  the  thin  layer  pro¬ 
perties  would  essentially  be  averaged  in  (on  a  thickness  weighted  basis), 
with  the  properties  of  the  surrounding  material  in  the  interpreted  model 
obtained  from  the  refraction  survey.  Although  this  appears  to  be  the  case 
in  certain  instances,  in  other  cases  the  presence  of  a  thin  layer  may  induce 
topography  or  fault-like  features  in  the  interpreted  subsurface  model  even 
though  the  true  model  consists  entirely  of  plane  layered  horizontally  lying 
geologic  structures. 

This  study  has  been  specific  rather  than  general  in  that  we  have  focused 
upon  a  subsurface  model  appropriate  for  the  alluvium-caliche  situation. 
Additional  work  will  be  required  to  develop  more  general  rules  regarding 
wavelength-layer  thickness  ratios  for  which  thin  layer  effects  may  occur. 

It  is  likely,  however,  that  the  results  obtained  apply  to  refraction  surveys 
conducted  in  any  area  where  the  wavelengths  are  large  compared  to  the  layer 
thicknesses.  The  presence  of  time  gaps  in  the  record  sections  (similar  to 
those  shown  in  Figure  1)  or  the  presence  of  extremely  emergent  first  arrivals 
(similar  to  those  seen  in  the  synthetic  record  sections)  should  alert  the 
interpreter  that  existence  of  thin  layer  effects  is  possible.  In  these  cases, 
the  results  of  the  refraction  surveys  should  be  corroborated  by  borehole 
data  or  other  independent  supporting  information.  Some  of  the  results  sug¬ 
gest  that  use  of  radial  geophones  for  refraction  surveys  may  be  advantageous 
in  the  detection  of  thin  layer  first  arrivals. 

The  results  of  this  study  also  have  implications  for  the  interpretation 
of  ground  shock  data  from  high-explosive  simulation  tests.  Ground  shock 
analysts  should  keep  in  mind  that  unexpected  effects  may  occur  in  the  pre¬ 
sence  of  geologic  layers  whose  thickness  is  relatively  small  compared  to  the 
wavelengths  of  the  ground  motions  propagating  through  the  test  bed.  In  some 
instances,  thin  layer  effects  might  be  erroneously  attributed  to  relief 
effects  or  other  phenomena. 


Increasing  the  frequency  of  the  source  pulse  might  increase  the  detecta¬ 
bility  of  the  thin  layer  and  reduce  the  severity  of  the  thin  layer  effects. 
This  appears  not  to  have  been  the  case,  at  least  for  the  moderately  higher 
frequency  source  used,  since  arrival  time  data  for  the  high  frequency  source 
shown  in  Figures  17  through  23  are  in  most  cases  indistinguishable  from  the 
lower  frequency  data.  Increasing  the  source  frequency  considerably  alters 
the  appearance  of  the  waveforms  as  shown  in  the  comparison  presented  in 
Figure  30  for  the  vertical  case  and  Figure  31  for  the  radial  case. 

Reference  18  presents  a  discussion  of  the  pervasiveness  of  the  hidden 
layer  or  blind  zone  problem  for  shallow  refraction  surveys  when  standard 
first  arrival  interpretation  techniques  are  used.  As  discussed  in  Reference 
18,  hidden  layer  problems  can  cause  serious  errors  in  the  interpretation  of 
shallow  refraction  surveys  for  engineering  site  characterization.  Hidden 
layer  problems,  coupled  with  the  thin  layer  effects  presented  in  this  report, 
suggest  the  need  for  a  more  robust  interpretation  method  for  shallow  seismic 
refraction  surveys.  Users  of  the  seismic  refraction  technique  should  be 
aware  of  its  limitations  and,  whenever  possible,  should  rely  upon  multiple 
sources  of  geophysical  data  for  construction  of  subsurface  geologic  models. 
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